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Abstract 


This  paper  describes  the  travel-time  anomalies  and 
attenuation  losses  of  seismic  compress i onal  waves  generated 
by  a  series  of  underwater  explosions  in  the  eastern  United 
States.  The  efficie.it  tamping  of  the  shots  fired  in  water 
provided  a  seismic  source  that  could  be  detected  at  much 
larger  ranges  than  could  be  accomplished  by  equivalent 
sized  shots  fired  underground. 

A  number  of  mobile  field  recording  stations  equippe 
with  three-component  matched  short  period  seismometers  and 
magnetic  tape  recorders  were  used  to  record  273  individual 
shots.  A  total  of  1,295  recordings  were  obtained  along  a 
reversed  profile  extending  from  International  Falls, 
Minnesota  to  the  Atlantic  Coast.  The  underwater  shots  were 
fired  in  Lake  Superior  and  the  Atlantic  Ocean.  Precise  tra¬ 
vel  times  were  obtained  by  recording  radio  time  signals  at 
all  of  the  recording  stations. 

An  analysis  of  the  travel-times  of  the  seismic  waves 
disclosed  that  the  earth's  crust  varies  in  thickness  from 
28.1  km  near  the  Atlantic  Coast  In  North  Carolina  to  50.4 
km  near  the  Keweenaw  Peninsula  In  upper  Michigan.  The  crust 
In  North  Carolin.  was  found  to  be  comprised  of  a  single  lay¬ 
er  with  a  compressional  wave  velocity  of  6.0  km/sec.  A  two- 
layer  crust  with  compressional  wave  velocities  of  6.2  and 
7.7  km/sec  was  disclosed  in  upper  Michigan.  Travel  time 
residuals  across  the  Appalachian  Mts.  would  indicate  a 
mountain  root  system  similar  to  that  found  under  the  Rocky 
Mts . 

The  Lake  Superior  area  was  found  to  be  more  efficient 
In  the  coupling  of  energy  1  ito  seismic  waves  by  the  under¬ 
water  shots  than  in  the  Atlantic  Ocean.  This  is  believed 
du:  mainly  to  the  bubble  pulse  phenomenon  anci  the  reinforce¬ 
ment  caused  by  reverberation  of  sound  waves;  between  the  lake 
bottom  and  the  water's  surface. 
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A  new  technique  of  displaying  seismic  attenuation  data 
In  three  dimensional  space  was  developed  that  permits  pre¬ 
sentation  of  large  quantities  of  data  In  a  concise  manner, 
hence  aiding  In  the  Interpretation  of  the  data.  Consistent 
anomalies  In  the  spectrum  of  the  seismic  waves  were  found 
that  could  be  correlated  with  various  parameters  in  the 
source  *,eg1on  and  with  propagation  effects. 

The  detailed  spectral  studies  that  were  made  as  a 
function  of  distance  from  the  source  disclosed  that  the 
attenuation  of  the  seismic  waves  could  be  expressed  by  the 
equation 

A  *  A  X'VafX 
o 

where  A  *  amplitude,  X  »  distance,  n  -  geometrical  spreading 
factor,  a  *  absorption  coefficient  and  f  =  frequency.  The 
crustal  attenuation  data  In  the  Lake  Superior  area  were  found 
to  have  a  value  of  n  *  1.5  and  a  »  2-3  x  10"^.  Upper  mantle 
compresslonal  wave  attenuation  data  gave  a  value  of  n  ■  1  and 
a  *  2  x  10“4. 
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Introduction 


The  attenuation  or  decrease  in  amplitude  of  seismic 
waves  with  increasing  distance  from  the  source  has  been  a 
problem  of  considerable  interest  to  seismologists  and  other 
scientists  for  many  years.  Predicting  damage  due  to  quarry 
shots,  areas  affected  by  earthquakes  and  the  general  problem 
of  detecting  and  locating  underground  nuclear  detonations 
are  just  a  few  areas  where  this  subject  has  specific  appli¬ 
cations.  Many  factors  can  contribute  to  the  seismic  loss 
of  energy  with  distance.  An  analysis  of  the  precise  travel 
times  of  seismic  waves  allows  seismologists  to  compute  the 
structure  of  the  earth's  interior  and  to  determine  the  seis¬ 
mic  velocities  of  the  structural  components.  Information 
concerning  the  latter  provides  insight  into  the  composition 
of  the  materials  that  comprise  these  features.  For  a  com¬ 
plex  problem  of  this  nature  sound  theoretical  hypotheses  and 
detailed  field  measurements  are  necessary  before  one  can 
ascertain  the  significant  contributing  factors.  It  is  the 
purpose  of  this  study  to  present  the  results  of  an  extensive 
investigation  on  the  propagation  of  seismic  waves  in  the 
eastern  United  States,  to  show  how  these  waves  are  attenuat¬ 
ed  and  to  determine  the  major  crustal  and  upper  mantle 
structure . 

Discussion  of  Wave  Types 

Factors  that  contribute  to  the  attenuation  of  eeismic 
waves  are  listed  below: 
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1.  absorption  ( 1 . e . ,  conversion  of  mechanical 
energy  into  heat) 

2.  scattering 

3.  geometrical  spreading 

4.  reflection  and  refraction  (including  conversion 
of  wave  types) 

Absorption  of  seismic  waves  s..ould  theoretically  be  of  a 
solid  friction  or  viscous  nature  or  a  combination  of  both. 

If  solid  friction  is  the  significant  factor,  an  exponential 
amplitude  decay  proportional  to  the  first  power  of  the  fre- 
quency  (A  -  e  where  a  *  absorption  coefficient  f  = 
frequency,  and  x  =  distance)  would  be  expected  (Howell  1963). 
Viscous  type  friction  on  the  other  hand  would  result  in  an 
exponential  decay  proportional  to  the  second  power  of  the 
frequency.  Born  [1941]  has  shown  fjr  low  frequency  data 
(below  150  cps)  that  viscous  absorption  is  unimportant  and 
if  scattering  occurs  from  small  inhomogeneities  less  than  a 
wavelength  In  dimension  that  the  loss  Is  of  the  same  nature 

-Otf  *  X 

as  viscous  absorption,  e.g.,  A  =  e  .  Selective  scatter¬ 

ing  due  to  reflections  and  refractions  on  prominent  Irregu¬ 
larities  would  have  more  pronounced  effects  on  wavelengths 
of  the  same  order  as  the  Irregularities.  Multiple  wave 
paths  may  cause  constructive  or  destructive  interference  re¬ 
sulting  in  scatter  of  the  attenuation  data. 

For  body  waves  (compressl onal  and  sh^ar)  geometrical 
spreading  would  cause  a  decrease  In  amplitude  that  would  be 
proportional  to  the  Inverse  first  power  of  the  distance 
(A  »  1 / x )  between  source  and  receiver.  Surface  waves 


3 


(Rayleigh  and  Love)  attenuate  like  the  inverse  one-half 
power  of  the  distance.  Head  v,ave  amplitudes  of  compres- 
sional  waves  (see  Heelan  [1953],  Werth  [1962]  or  Donato 
[1964])  should  fall  off  like  (x")‘2  where  x'  =  (xL3)1/4 
with  x  equal  to  the  horizontal  distance  from  the  source 
to  the  recording  station  and  L  the  distance  traveled  in 
the  refractor.  Hence  for  large  distances  head  waves 
should  attenuate  approximately  like  x"2.  The  last  major 
factor  contributing  to  the  attenuation  of  seismic  waves  is 
the  energy  loss  that  a  wave  encounters  upon  reflection  and 
refraction  (including  conversion  of  wave  type)  at  a  ve¬ 
locity  discontinuity.  This  factor  has  been  described  in 
some  detail  by  the  classical  works  of  Knott  [1899], 

Macelwane  [1936],  Zoeppritz  [1919],  and  others. 

This  study  will  be  restricted  primarily  to  the  atten¬ 
uation  of  body  waves  of  the  compressi onal -di 1 atati onal  type. 
Travel  time  data  will  be  included  to  provide  as  much  corre¬ 
lation  with  geological  parameters  as  possible.  The  follow¬ 
ing  section  contains  a  description  of  the  field  data  collect¬ 
ed  over  a  five  year  t-'me  period  that  was  utilized  in  this 
investigation. 


Field  Measurement  Program 

A  reversed  long-range  seismic  refraction  profile  was 
conducted  along  a  northwest-southeast  line  approximately 
2000  km  long  in  the  eastern  United  States.  The  line  extended 
from  International  Falls,  Minnesota  to  the  North  Carolina 


Atlantic  Coast.  This  exercise  comprised  a  series  of  separate 
experiments  with  many  participants.  The  data  collected  by 
University  of  Michigan  field  crews  is  unique  in  that  it  is 
one  of  the  longest  reversed  profiles  with  close  station  con¬ 
trol  in  existence,  and  It  is  also  the  only  reversed  profile 
over  the  Appalachian  Mountains.  The  following  series  of 
shots  comprise  this  study:  (see  Table  I). 

Details  of  these  experiments  are  Included  in  Appendix 
A.  The  University  of  Michigan  operated  from  two  to  nine 
portable  seismograph  stations  to  record  these  shots.  Each 
station  was  equipped  with  FM  magnetic  tape  recorders,  match¬ 
ed  short-period  three-component  seismometers  and  short  wave 
radios  to  record  radio  time  signals.  Details  and  specifica¬ 
tions  of  this  equipment  are  discussed  In  Appendix  B.  The 
general  technique  employed  to  obtain  the  detailed  close  sta¬ 
tion  control  was  to  assign  mobile  stations  to  a  given  seg¬ 
ment  along  the  main  refraction  profile.  The  shots  were 
usually  fired  at  night  ranging  over  a  time  period  from 
several  days  to  a  month.  The  recording  stations  were  moved 
each  day  so  as  to  cover  as  much  ground  within  their  segment 
as  possible.  Sites  were  located  on  bedrock  wherever  avail¬ 
able  and  in  areas  to  avoid  cultural  activity.  At  the  larger 
distances  with  low-signal  levels  the  site  selection  was  very 
critical  since  poor  signal  *  .o-nolse  levels  made  it  very  dif¬ 
ficult  to  detect  the  first  compressiunal  wave  arrivals. 


Table  I.  Explosion  Seismic  Experiments 
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Travel  Times  and  Crustal  Structure 


The  shot  point  and  recording  station  locations  are 
shown  in  Figure  1  where  data  reported  in  this  paper  were  ob¬ 
tained.  The  series  of  shots  fired  in  the  Atlantic  Ocean  in¬ 
clude  the  CHASE  series,  the  1962  North  Carolina  Experiment, 
and  the  ECOOE  project.  The  Lake  Superior  shots  mentioned 
in  the  previous  section  are  also  shown  in  this  figure. 

Reduced  travel  time  graphs  were  prepared  for  the  first 
compressional  wave  arrivals  recorded  at  the  stations  shown 
it.  this  map.  These  are  presented  in  Figures  2  through  12. 
Graphs  of  this  type  allow  one  to  examine  more  carefully  the 
deviations  or  anomalies  in  travel  times  than  would  be  pos¬ 
sible  to  see  in  the  standard  time  vs  distance  graphs. 

Figure  2  is  a  composite  made  up  of  first  compressional 
wave  arrivals  recorded  at  Jacksonville,  Potters  Hill,  and 
Dudley,  North  Carolina  for  the  shots  fired  during  the  1962 
North  Carolina  experiment.  The  University  of  Michigan  sta¬ 
tion  locations  and  first  arrival  times  may  be  found  in  a 
paper  by  Meyers,  et  al .  [1965].  Using  a  crustal  velocity 
of  6.0  km/sec  for  the  reduced  travel  time  curves  it  can  be 
seen  that  velocities  of  5.8  and  7.9  km/sec  for  the  crustal 
layer  and  upper  mantle  velocity  gives  a  reasonable  fit  to 
the  data.  Figures  3  and  4  show  the  reduced  travel  time  re¬ 
siduals  for  the  individual  stations  of  the  ECOOE  Southern 
Profile  shots.  The  station  locations  and  first  arrival 
times  for  the  150,  250,  and  400  km  University  of  Michigan 
stations  on  the  Southern  Profile  and  the  150,  300,  and  325 
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FIGURE  2.  TRAVEL-TIME  RESIDUALS:  1962  NORTH  CAROLINA  EXPERIMENT 
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FIGURE  3.  TRAVEL-TIME  RESIDUALS:  SOUTHERN  PROFILE 
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FIGURE  4.  TRAVEL-TIME  RESIDUALS:  SOUTHERN  PROFILE,  400- , 
450- ,  550- ,  600- ,  650- ,  AND  750- kn  SITES. 


km  stations  on  the  Northern  Profiles  of  the  ECOOE  series 
may  be  found  in  a  paper  by  Hales,  et  al  .  [1967].  This  in¬ 
formation  was  supplied  by  the  University  as  a  part  of  theiir 
contribution  to  the  cooperative  experiment  conducted  by  the 
North  American  Explosion  Seismology  Group.  Figure  5  is  a 
composite  of  the  data  for  Figures  2  through  4  and  include 
the  CHASE  III  and  VII  shots.  An  upper  mantle  velocity  of 
8.2  km/sec  was  used  to  compute  the  data  shown  in  the  latter 
figures . 

The  data  shown  in  Figure  5  disclose  three  branches  in 
the  travel  time  curve.  These  may  be  fit  reasonably  well  to 
the  following: 

T-j  =  a/2.0  +  0  (0-30  km)  h^  =  0.49  km  (assumed)  (1) 

T2  ■  a/6.0  +  0.48  (30-150  km)  h?  =  28.69  (2) 

T3  =  a/8.2  +  7.0  (150-1000  km)  (3) 

where  T  -  travel  time  of  first  compressi onal  wave  arrival 

A  -  distance  between  shot  point  and  recording  station 
h  =  layer  thickness. 

The  velocity  for  the  first  branch  or  upper  layer  is 
assumed  since  we  have  no  data  at  close  ranges.  A  velocity 
of  2.0  km/sec  is  quite  reasonable  for  a  surface  layer.  An 
error  in  this  assumption  would  produce  no  serious  discrepan¬ 
cies  for  the  total  thickness  of  the  crust  since  the  inter¬ 
cept  of  the  second  branch  (C.43  sec)  is  small.  These  data 
agree  quite  well  with  those  reported  by  Hales  which  will  be 
discussed  later. 
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Considerable  scatter  is  evident  on  the  upper  mantle 
(8.2  km/sec)  branch  of  the  travel  time  curve  shown  in  Figure 
5.  However,  when  observing  the  individual  station  residuals 
in  Figures  3  and  4  a  systematic  deviation  can  be  seen.  The 
150,  250,  and  350  km  sites  have  a  residual  of  approximately 
7  to  7.5  seconds.  The  latter  two  sites  indicate  a  slightly 
higher  Pn  velocity.  The  400  km  site  which  is  located  in 
the  Piedmont  Belt  shows  a  significantly  lower  residual  of 
approximately  6  seconds.  This  could  be  accounted  for  by  a 
thinner  crustal  section  related  to  the  Piedmont  Triassic 
Belt.  The  western  edge  of  the  Appalachian  root  lies  approx¬ 
imately  halfway  between  the  450  and  550  km  sites.  An  in¬ 
crease  in  residuals  from  7.5  to  8.5  seconds  seen  at  these 


two  sites  could  be  caused  by  the  effect  of  the  root.  This 
would  also  eff°ct  the  next  two  sites  (600  and  650  km)  which 
would  fall  in  the  shadow  of  the  root.  The  stations  at 


greater  distances  should  not  show  any  anomalous  effects  of 
the  root  since  the  ray  paths  would  be  below  the  projection 
of  the  root  system  which  is  reported  by  James  [1966]  to 
reach  60  km.  As  can  be  seen  in  Figure  5  the  residuals  be¬ 
yond  700  km  average  about  7.5  seconds. 

Similar  travel  time  residuals  for  the  EC00E  Northern 
Profile  are  shown  in  Figures  6  through  8.  Unfortunately  the 
shot  point  location  accuracy  of  the  shots  fired  off  the  con¬ 
tinental  shelf  are  uncertain.  The  shots  are  normally  located 
by  using  water  wave  arrivals  recorded  at  underwater  hydro¬ 


phones.  For  the  Northern  Profile  shots  the  hydrophone  buoys 
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could  n.t  be  placed  in  proper  locations  so  that  it  was 
necessary  to  rely  on  a  combination  of  hydrophone  arrays 
and  LORAN  C.  In  analyzing  the  data  afterwards  the  South¬ 
west  Graduate  Research  Center  found  that  rather  large  dis¬ 
crepancies  existed  between  locations  determined  by  the  two 
techniques  (LORAN  C  had  errors  of  up  to  10%  due  in  part  to 
the  fact  that  part  of  the  transmission  path  was  over  land). 
This  uncertainty  in  shot  location  can  be  seen  in  the  range 
of  travel  time  residuals  for  shots  20  through  37.  Also 
contributing  to  the  scatter  in  residuals  is  the  fact  that 
these  shots  (20  through  37)  w«;re  fired  in  deep  water  (up 
to  12,240  ft)  at  a  depth  of  300  feet  while  the  other  shots, 
located  on  the  continental  shelf,  were  fired  on  the  ocean 
bottom.  All  of  the  data  should  be  corrected  to  a  common 
datum  taking  into  account  the  difference  in  velocities  be¬ 
tween  salt  water  and  the  sediments  on  the  continental  shelf. 
Since  this  particular  problem  is  not  germane  to  the  main 
topic  of  this  report  and  since  the  correct  velocities  to  be 
used  in  this  adjustment  are  not  known  precisely,  no  attempt 
will  be  made  to  correct  the  travel  time  residuals  for  the 
shots  fired  off  the  shelf.  The  residuals  on  the  Pn  branch 
of  the  travel  time  curve  for  the  150,  300,  325,  350,  and 
375  km  sites  average  between  7.5  to  8  seconds  which  would 
be  indicative  of  a  slightly  thicker  crust  than  observed  on 
the  Southern  Profile.  The  western  edge  of  the  Appalachian 
root  is  near  the  350  km  site.  The  stations  at  400,  425, 
and  450  km  show  an  increase  in  this  residual  of  about  1  sec 
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which  verifi.s  the  shadow  effect  of  the  root  which  was  ob¬ 
served  on  the  Southern  Profile.  No  evidence  of  a  thinner 
crustal  section  which  was  observed  in  the  Piedmont  section 
of  the  Southern  Profile  was  found  on  this  line.  However, 
an  examination  of  the  geology  of  this  region  shows  that  the 
width  of  the  Appalachian  Mountain  system  is  considerably 
narrower  in  this  area  and  that  the  Piedmont  section  is  quite 
thin.  Hence,  one  could  not  expect  a  similar  seismic  anoma¬ 
ly  on  this  prof i 1 e . 

Figure  8  is  a  composite  of  all  sites.  The  residuals 
for  shots  22  through  37  have  been  omitted  on  this  graph. 

The  larger  average  positive  anomaly  between  400  and  500  km 
Is  quite  apparent.  However,  the  large  anomaly  at  880  to 
970  km  which  was  recorded  at  the  Ann  Arbor  Botanical  Gardens 
Station  is  rather  surprising.  A  station  anomaly  of  this 
magnitude  was  not  substantiated  by  the  CHASE  IV  or  VII  shots. 
However,  the  CHASE  III  shot  recorded  at  the  Botanical  Gar¬ 
dens  and  at  a  station  160  km  from  Ann  Arbor  seem  to  sub¬ 
stantiate  it.  Two  possible  explanations  can  be  offered. 

One  Is  that  the  events  identified  as  the  Pn  arrival  were  in¬ 
correctly  picked  and  what  was  identified  as  Pn  was  actually 
a  later  larger  amplitude  P  arrival.  The  other  possibility 
could  be  a  shot  point{s)  anomaly.  The  latter  does  not 
appear  very  feasible  since  CHASE  III  and  IV  were  fired  in 
the  same  general  area. 

The  data  on  the  Northern  Profile  may  be  divided  into 
three  groups  and  can  be  expressed  mathematically  as: 
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h1  =  1.68 

(4) 

=  A/ 6.0  +  1.63 

h2  =  25.84 

(5) 

=  A/8.2  +  7.5 

(6) 

All  of  th?  participants  including  the  University  of 
Michigan  submitted  their  first  arrival  times  on  the  ECOOE 
project  to  the  Southwest  Graduate  Research  Center.  Hales 
[1967]  has  analyzed  these  data  and  reports  the  following. 
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Southern  Profile  (normal  to  Coast) 


T  =  A/1 .70  +  0 

hi 

=  0.49 

(7) 

T  =  A/6.03  +  0.55 

h2 

=  30.38 

(8) 

T  =  a/8.13  +  7.32 

(9) 

Northern  Profile 

(normal 

to  Coast) 

T  =  a/2.10  +  0 

hi 

=  1.63 

(10) 

T  =  A/5.78  +  1 .45 

h2 

=  8.31 

(ID 

T  =  a/6.34  +  2.65 

h3 

=  16.32 

(12) 

T  =  A/7.97  +  6.60 

(13) 

Figure  9  shows  the  combined  residuals  for  the  five- 
and  ten-ton  shots  of  the  ECOOE  series  and  all  of  the  CHASE 
shots.  The  data  shown  in  this  figure  and  in  Figure  5  yield 
the  following. 
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THE  ECOOE  SERIES  AND  CHASE  III,  IV,  AND  VII. 
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A  least-squares  analyses  was  made  for  the  travel-time 
residuals  shown  for  each  station  in  Figures  3,  4,  6,  and  7. 
The  deep  water  shots  were  omitted  because  of  the  large  vari¬ 
ations  in  residuals.  The  resultant  velocities  and  inter¬ 
cepts  are  shown  in  Table  II.  The  data  for  those  stations 
which  have  only  a  few  recordings  or  where  the  distance  range 
is  limited  should  be  used  with  caution. 


Table 

II . 

ECOOE  Compressi onal 

Veloci ties 

Northern  Profile 

Site 

Velocity  Intercept 

( km/sec) 

(sec) 

1  50 

I 

6.8 

3.89 

150 

II 

8.2 

7.48 

150 

III 

8.0 

6.95 

300 

8.0 

6.69 

325 

8.1 

7.42 

350 

3.7 

10.29 

375 

8.0 

6.56 

40C 

7.7 

5.97 

425 

8.2 

8.82 

450 

8.4 

10.52 

600 

Insufficient 

data 

925 

8.5 

13.99 

Southern  Profile 

150 

I 

6.0 

0.57 

150 

II 

8.0 

6.81 

250 

8.5 

8.36 

350 

8.7 

9.56 

400 

8.7 

9.20 

450 

8.2 

7.30 

550 

8.5 

10.80 

600 

9.0 

14.69 

650 

9.1 

16.71 

750 

Insufficient 

data 

The  reversed  profile  for  the  series  of  shots  fired  in 
the  Atlantic  Ocean  is  made  up  of  a  composite  of  shots  fired 
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in  Lake  Superior.  These  shots  include  the  1963  Lake 
Superior  Experiment,  the  follow-up  sequence  of  shots  in 
1964  by  The  University  of  Wisconsin,  the  October  1964  series 
of  shots  conducted  by  the  U.  S.  Geological  Survey  end  the 
1966  Project  Early  Rise  shots. 

The  1963  Lake  Superior  Experiment  consisted  of  a  series 
of  38  lake  bottom  shots  fired  every  10-15  km  along  the 
length  of  Lake  Superior.  The  University  of  Michigan's  par¬ 
ticipation  in  this  experiment  is  described  by  Willis,  et  al  . 
[1964],  Steinhart  [1964],  Willis  and  DeNoyer  [1966]  and  in 
other  papers  by  the  author.  Five  major  profiles  were  ob¬ 
tained  from  this  experiment  only  one  of  which  will  be  utiliz¬ 
ed  in  this  report.  A  profile  extending  southeastward  from 
Houghton  on  the  Keweenaw  Peninsula  in  upper  Michigan  was 
obtained  by  using  three  mobile  stations  to  occupy  13  record¬ 
ing  sites.  Shot  points  located  opposite  Houghton  were  re¬ 
corded  at  distances  ranging  from  84  to  258  km.  The  reduced 
travel  time  data  from  this  profile  are  shown  in  Figure  10. 
Some  scatter  is  evident  that  appears  related  to  individual 
station  corrections.  A  crustal  velocity  of  6.4  km/sec  was 
obtained.  A  lower  crustal  velocity  of  7.5  km/sec  was  also 
indicated.  However,  the  scatter  of  the  data  made  this  value 
somewhat  doubtful . 

The  July  1964  Lake  Superior  Experiment  did  not  contri¬ 
bute  significantly  to  this  study  since  the  recording  stations 
and  shots  were  located  to  supplement  the  Eastern  and  Western 
Profiles  on  the  1963  Lake  Superior  Experiment.  The  shots 
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FIGURE  10.  TRAVEL-TIME  RESIDUALS  FOR  1963  LAKE  SUPERIOR  SOUIH  LINE 
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of  the  October  1964  series  were  utilized  to  augment  and 
extend  the  1963  South  Line  to  Ann  Arbor,  Michigan  (see 
Willis  and  Tanis  [1965]).  These  shots  consisted  of  1,  2, 

3,  6,  and  10  tons  and  were  fired  at  several  shot  point  lo¬ 
cations.  At  least  one  shot  of  each  size  was  detonated  near 
the  location  of  snot  point  45  of  the  1963  series.  The  re¬ 
duced  travel  time  graph  for  this  series  is  shown  in  Figure 
11.  A  crustal  velocity  of  6.17  km/sec  and  a  Pn  velocity  of 
8.29  km/sec  gave  a  reasonable  fit  to  these  data. 

When  plans  were  made  for  participating  in  Project  Early 
Rise  it  was  decided  to  duplicate  or  reoccupy  the  stations 
used  on  the  1963  Lake  Superior  South  Line  and  the  1964 
October  profile  as  well  <s  extending  the  line  southeastward 
to  the  Atlantic.  Several  reasons  prompted  this  decision. 

Better  control  was  needed  to  evaluate  some  early  arrivals 
detected  nea.'  the  intersection  of  the  crustal  and  mantle 
branches  of  the  travel  time  curve;  the  variation  in  charge 
size  of  the  1964  shots  made  difficult  the  interpretation  of 
the  attenuation  data;  and  finer  station  spacing  was  desirable 
-■'cross  the  Michigan  Basin. 

All  of  the  Early  Rise  shots  were  fired  in  the  vicinity 
of  shot  point  51  of  the  1963  series.  The  reduced  travel 
time  data  for  the  vertical  component  of  the  fi^st  compressional 
wave  arrivals  for  the  Early  Rise  shots  recorded  along  the 
main  northwest-southeast  line  are  shown  in  Figure  12.  An 
upper  mantle  velocity  of  8.2  km/sec  was  used  for  these  com¬ 
putations.  A  crustal  velocity  of  6.2  km  veri.  s  the  results 
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FIGURE  11.  TRAVEL-TIME  RESIDUALS  FOR  OCTOBER  1964  U\KE  SUPERIOR 
EXPERIMENT. 


FIGURE  12.  TRAVEL-TIME  RESIDUALS  -  OR  EARLY  RISE. 
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shown  in  Figure  11.  A  clear-cut  lower  crustal  layer  with 
a  velocity  of  7.7  km/sec  was  obtained  verifying  the  inter¬ 
mediate  layer  that  was  suspected  from  the  data  shown  in 
Figures  1 0  and  1 1 . 

Considerable  scatter  can  be  seen  in  the  Pn  arrivals 
beyond  700  km.  This  is  due  primarily  to  the  high  seismic 
background  r.oise  caused  by  increased  cultural  activity 
across  southeastern  Michigan,  Ohio,  and  West  Virginia.  The 
difficulty  ir.  determining  the  first  arrivals  along  this 
profile  can  be  appreciated  by  examining  Figure  16c,  page  97, 
of  a  report  by  Warren,  et  al.  [1967]  which  shows  a  composite 
record  section  of  this  line.  The  details  at  the  scale  used 
are  difficult  to  see.  However,  the  increased  noise  beyond 
750  km  Is  quite  evident.  One  feature  particularly  inter¬ 
esting  Is  the  small  amplitude  first  arrivals  at  stations 
located  at  160  to  275  km.  This  is  illustrated  better  in 
Figure  13.  At  these  distances  the  record  gains  in  the  fiel‘ 
were  deliberately  turned  up  on  the  high  gain  channels  so 
that  the  record  would  be  overrecorded.  This  provided  better 
record  gains  for  the  low  amplitude  first  arrivals.  This 
figure  shows  vhe  first  arrivals  on  the  vertical  component 
with  the  overrecorded  later  phases  omitted.  The  records  are 
adjusted  so  that  arrivals  with  a  velocity  of  7.7  km/sec 
should  be  aligned  vertically  on  this  section.  Fhe  top  three 
traces  are  on  the  6.2  km/sec  branch  of  the  travel  time  curve. 
If  the  record  gains  used  in  the  field  were  set  so  the  later 
phases  would  not  be  overrecorded,  this  early  arrival  would 
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be  very  difficult  if  not  impossible  to  detect.  There  is  a 
second  recording  of  another  shot  at  each  of  the  stations 
shown  here.  These  second  shots  make  almost  perfect  over¬ 
lays  with  those  shown  in  this  figure  so  there  can  be  no 
doubt  about  the  validity  of  the  apparent  intermediate  lay¬ 
er.  This  layer  correlates  with  the  Conrad  discontinuity  re¬ 
ported  in  Europe. 

Figure  14  is  a  composite  tracing  of  the  vertical 
component  first  arrivals  for  all  of  the  Early  Rise  shots  re¬ 
corded  at  the  University  of  Michigan  Botanical  Gardens  in¬ 
termediate-depth  well  station.  The  record  gain  used  (in  db) 
for  each  shot  is  shown  along  the  right  hand  margin  of  the 
figure.  The  repeatability  of  the  wave  form  is  quite  good. 
Varying  background  noise  prior  to  the  first  arrivals  would 
make  the  identification  of  the  first  arrival  difficult  in 
some  cases  where  only  one  or  two  recordings  were  available. 

For  this  reason  at  some  of  the  more  distant  locations  four 
or  more  shots  were  recorded  in  an  effort  to  increase  the 
accuracy  of  the  arrival  times  and  ground  amplitudes. 

The  evidence  for  an  Appalachian  Mountain  root  system 
disclosed  on  the  ECOOE  Southern  Profile  is  not  as  convincing 
on  this  profile.  The  mountains  are  at  a  distance  of  1175  to 
1360  km  from  the  shot  poult.  The  residuals  for  those  stations 
located  on  the  western  half  of  the  mountains  are  smaller  than 
those  on  the  center  and  eastern  portion  and  for  those  sta¬ 
tions  located  just  to  the  east  of  the  mountains.  However, 
the  variation  in  the  data  closer  to  the  shot  point  would 
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FIGURE  14.  COMPOSITE  SEISMOGRAM  OF  EARLY  RISE  SHOTS 

RECORDED  AT  ANN  ARBOR. 
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cause  one  to  think  this  anomaly  could  also  be  due  to  just 
scatter  in  the  data.  For  the  shots  recorded  at  distances 
of  1500  to  1650  km,  the  first  arrival  could  not  be  read 
accurately.  At  our  most  distant  station  the  signal  ampli¬ 
tude  was  larger  even  though  the  background  noise  was  high. 

The  vertical  component  recordings  for  four  shots  recorded 
at  this  furthest  station  are  shown  in  Figure  15.  These  four 
individual  shots  were  dubbed  onto  one  magnetic  tape  taking 
into  account  shot  zero  times  and  slight  differences  in  shot 
point  locations.  The  times  were  aligned  so  that  each  shot 
was  within  ±2  milliseconds  at  the  time  of  the  first  arrival. 
Th<»  maximum  deviation  over  the  entire  record  is  ±5  milli¬ 
seconds.  The  vertical  components  shown  here  were  filtered 
on  playback  with  a  1  to  4  cps  passband  which  greatly  improved 
the  signal -to-noise  ratio.  The  summed  vertical,  longitudi¬ 
nal  and  transverse  components  for  the  four  shots  are  also 
shown  in  this  figure.  At  this  distance  it  can  be  seen  that 
most  of  the  energy  of  the  first  comprcssional  wave  is  re¬ 
corded  on  the  vertical  component.  The  summing  technique 
(type  of  array  processing)  and  bandpass  filtering  gives  a 
12  db  improvement  in  signal  to  noise  for  the  first  arrival. 
Such  techniques  allow  a  more  accurate  determination  of  the 
first  arrival  times  but  It  Is  quite  laborious  to  make  the 
composite  magnetic  tape  dubs.  This  process  was  used  only 
where  control  points  were  needed  for  critical  interpretations 
and  where  four  or  more  shots  were  recorded  at  the  same  sta¬ 


tion  . 
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FIGURE  15.  EARLY  RISE  SHOTS  lft-21  °ECOROED  AT  POTTER'S  HILL,  NORTH  CAROLINA 
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The  combined  travel  time  data  of  Figures  10-12  yield 
the  following: 


T]  =  A/4.8*  +  0 

h]  =  2.9 

(18) 

T2  =  A/6.2  +  1.0 

h2  =  22.5 

(19) 

T  3  =  A/7.7  +  5.9 

h3  =  25.0 

(20) 

T^  =  A/8.3  +  8.5 

h4  =  21.6 

(21) 

T,  =  A/8.4  +  9.3 

(22) 

*From  Smith,  et  al  .  (1  966) 

A  comparison  of  the  crustal  and  upper  mantle  structure 
at  both  ends  of  this  long  range  refraction  profile  shows  con¬ 
siderable  differences.  The  total  crustal  thickness  at  the 
continental  margin  is  28.1  km  compared  to  50.4  km  at  Lake 
Superior.  Two  major  crustal  layers  with  velocities  of  6.2 
and  7.7  km/sec  are  clearly  evident  at  the  latter  area  while 
only  a  single  layer  with  a  velocity  6.0  km/seo  was  found  at 
the  continental  margin.  A  thin  layer  of  low  velocity  sedi¬ 
ments  at  the  surface  is  assumed  in  both  areas.  A  slightly 
higher  upper  mantle  velocity  (8.3  vs  8.2  km/sec)  was  found 
in  the  Lake  Superior  area  while  the  higher  velocity  layer 
in  the  mantle  was  found  at  a  shallower  depth  (72  km  vs  99  km 
at  the  continental  margin).  However,  the  poor  signal-to- 
noise  ratios  and  the  scatter  of  the  data  make  the  validity 
of  the  latter  somewhat  questionable. 

If  the  data  shown  in  Figure  11  were  replotted  in  Figure 
12  using  the  same  upper  mantle  velocity  of  8.2  km/sec  for 
computing  the  residuals,  the  Pn  branch  of  the  travel  time 
residuals  for  the  October  1964  shots  would  average  about  1 


second  later  than  those  of  the  Early  Rise  shots.  The  6.2 
km/sec  branch  agrees  quite  well  while  the  7.7  km/sec  branch 
is  approximately  0.3  seconds  late.  This  could  only  be  ex¬ 
plained  by  a  difference  in  crustal  structure  between  shot 
points  45  and  51.  Figure  22  in  Smith,  et  al.  (1966)  shows 
that  shot  point  45  is  located  on  the  flank  of  an  anomalous 
rise  in  the  surface  of  the  upper  mantle  and  that  a  trough 
or  depression  in  the  mantle  is  apparent  between  these  two 
shot  points.  The  observed  discrepancy  of  one  second  would 
indicate  a  difference  in  crustal  thickness  of  the  order  of 
6  or  7  km.  The  mantle  depression  reported  by  Smith,  et  al  . 
(1966)  is  not  quite  that  large.  However,  a  ray  traversing 
this  feature  obliquely  would  not  require  the  depression  to 
be  as  deep  as  6  km  to  cause  the  one  second  travel  time  dis¬ 
crepancy  . 

Many  other  authors  (Smith,  et  al.  [1966],  Berry  and 
West  [1966],  Mereu  [1966],  Cohen  and  Meyer  [1966],  and  others) 
report  a  single  layer  crust  with  a  high  velocity  (6.5  to  7.2 
km/sec).  Our  Copper  Harbor  East  and  West  profiles  for  the 
1963  Lake  Superior  Experiment  showed  only  one  crustal  layer 
with  a  velocity  of  6.60  to  6.71  km/sec.  However,  the  Knife 
River,  Minnesota  Profile  disclosed  two  crustal  layers  with 
velocities  of  6.19  and  7.01  km/sec.  Green  and  Hales  [1967] 
and  Roller  and  Jackson  [1966]  report  a  two-layer  crust  on 
their  southwestern  profiles  from  Lake  Superior  with  veloci¬ 
ties  of  6. 3/5. 5  and  6.97/7.2  km/sec  respectively.  Why  then 
should  our  profile  extending  southeastward  from  the  Keweenaw 
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Peninsula  show  the  two-layer  crust  while  our  other  two 
Keweenaw  Peninsula  profiles  indicate  only  one?  The  simplest 
explanation  lies  with  the  fact  that  the  onset  of  the  7.7 
km/sec  arriva1  has  a  very  small  amplitude  (a  factor  of  30 
lower  than  the  later  arrival  corresponding  to  6.2  layer). 

This  lower  amplitude  event  could  go  undetected  because  where 
only  one  record  gain  is  used,  that  gain  would  be  set  so  as 
not  to  overrecord  the  later  arrival.  We  used  double  or  tri¬ 
ple  gains  to  investigate  this  anomaly.  For  the  Early  Rise 
shots  we  used  close  station  separations  at  distances  of  140 
to  280  km  to  verify  the  suspected  two-layered  crust.  The 
fact  that  the  same  shot  point  was  used  each  time  also  great¬ 
ly  facilitated  this  study  since  the  small  amplitude  first 
arrivals  could  be  identified  by  overlaying  the  records  re¬ 
corded  at  the  same  stations.  Such  was  not  possible  for  the 
1963  series. 

The  total  thickness  of  the  crust  is  not  altered  signi¬ 
ficantly  (only  0.1%)  by  the  two-layered  crust  with  velocities 
of  6.2  and  7.7  km/sec  compared  to  a  single-layered  crust  of 
6.8. 


Attenuation  Studies 

The  detonation  of  large  high  explosive  shots  in  several 
hundred  feet  of  water  provides  a  very  efficient  coupling 
mechanism  for  seismic  waves  (see  Weston  [1960],  DeNoyer  and 
Willis  [1961]  and  Willis  [1963]).  The  increased  signal 
level  for  a  properly  tamped  water  shot  can  generate  ground 
displacements  that  are  a  factor  of  20  or  more  larger  than 


an  equivalent  sized  shot  fired  in  ruck.  Hence,  the  firing 
of  1  to  10  ton  shots  in  Lake  Superior  and  off  the  Atlantic 
Coast  has  provided  seismologists  with  a  source  that  can  be 
detected  to  ranges  on  the  order  of  2500  km  under  favorable 
condi ti ons . 

Earlier  research  (Willis,  1964)  had  shown  that  the 
attenuation  of  high  frequency  (above  5  cps)  seismic  energy 
at  distances  of  several  hundred  kilometers  was  significant¬ 
ly  lower  in  the  northeastern  United  States  than  in  the 
western  areas  of  the  country.  Hence  recording  equipment 
was  used  that  would  insure  recording  any  high  frequency  en¬ 
ergy  that  might  be  generated  by  these  underwater  shots. 

EC00E  Series 

The  observed  loss  in  amplitude  as  a  function  of  distance 
from  the  1962  North  Carolina  shots,  the  EC00E  series,  the 
CHASE  shots,  and  the  1963-64  Lake  Superior  shots  has  been 
reported  in  nine  previous  papers  by  the  author.  In  Figures 
15-18,  21,  24,  and  27  of  Willis  and  Jackson  [1966]  the  ver¬ 
tical  component  of  the  first  compressional  wave  arrival  for 
the  EC00E  shots  car  be  seen  to  fall  off  like  the  inverse  cube 
of  the  distance.  Figures  22  and  25  likewise  show  a  si,.  :lar 
fall  off  for  the  maximum  amplitude  P  wave  arrival.  The  maxi¬ 
mum  surface-wave  arrival  shown  in  Figures  23  and  26  can  be 
seen  to  fall  off  at  a  rate  between  the  inverse  square  and 
the  inverse  cube  of  the  distance. 


Pn  attenuation  curves  for  the  individual  5-ton  shots 
of  the  ECOOE  northern  profile  are  shown  in  Figure  16.  The 
major  fluctuation  in  amplitude  correlates  with  those  sta¬ 
tions  in  the  Appalachian  Mountains.  The  signal  levels  of 
the  vertical  component  of  first  compres s i ona 1  wave  (P)  ar¬ 
rival  and  the  maximum  surface  wave  (R}  arrival  are  shown 
in  Figure  17  for  the  ECOOE  10-ton  shots  (labeled  6  and  39) 
of  the  southern  profile  and  for  CHASE  III  and  IV  (labeled) 

3  and  4).  The  rate  cf  attenuation  is  approximately  propor¬ 
tional  to  the  inverse  cube  of  the  distance.  An  anomalous 
increase  in  signal  with  distance  was  observed  across  the 
Michigan  Basin.  It  is  also  interesting  to  note  chat  the 
maximum  surface  wave  amplitudes  for  the  10-ton  shots  (6R 
and  39R )  a>":  of  the  same  order  of  magnitude  at  many  distances 
as  the  compress i onal  wave  amplitudes  for  CHASE  III  and  IV. 

For  the  same  shot  the  maximum  amplitude  surface  waves  ranged 
10  to  20  db  higher  than  the  first  arrival  compressional  waves. 

October  1964  Lake  Superior  Shots 

Attenuation  studies  of  this  series  of  recordings  are 
complicated  by  the  fact  that  the  charge  size  varied  from  cne 
to  ten  tons  and  multiple  shot  points  were  used.  Comparing 
recordings  made  for  the  same  shot,  disclose  anomalou'  ampli¬ 
tude  variations  such  as  that  shown  in  Figure  18.  F.iown  here 
are  the  vertical  component  ground  particle  velocity  ampli¬ 
tude:  of  the  first  compressional  wave  arrival  for  five  sta¬ 
tions  ranging  in  distance  from  549  km  to  766  km.  The  two 
most  distant  stations  had  larger  peak  amplitudes  than  the 
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FIGURE  16.  SIGNAL-LEVEL  vs  DISTANCE  FOR  B-^ON  SHOTS  OF  ECOOE 
NORTHERN  PROFILE. 
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FIGURE  17.  SIGNAL- LEVEL  vs  DISTANCE  FOR  10-TON  SHOTS  ECOO! 
SOUTHERN  PROFILE  AND  CHASE  III  AND  IV. 
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FIGURE  18.  SPECTRAL  ANALYSIS  OF  10-TON  LAKE  SUPERIOR  SHOT  RECORDED 
AT  STATIONS  LOCATED  IN  MICHIGAN  BASIN. 


close  stations  and  the  center  static,,  had  the  lowest  ampli¬ 
tude.  Some  of  these  variations  are  probably  due  to  local 
geological  factors.  However,  all  of  the  stations  are  lo¬ 
cated  in  the  -  ame  type  of  geological  environment  (southeast 
flank  of  the  Michigan  Basin)  and  are  spaced  relatively  close 
together . 

In  order  to  determine  a  mathematical  expression  that 
would  best  fit  all  of  the  recorded  first  compress ional  wave 
data  from  this  experiment  a  digital  computer  program  was 
written  that  would  compute  the  essential  variables.  In  or¬ 
der  to  eliminate  shot  point  effects,  only  those  recordings 
from  shots  fired  at  shot  point  45  were  used. 

We  assume*  that 


where  A  =  particle  velocity  amplitude 
x  =  distance 

n  =  geometrical  spreading  factor 
a  =  absorption  coefficient 
f  =  frequency 
W  =  charge  size 

Solving  for  a,  we  have 


*See  Willis  [1960],  Willis  and  Wilson  [1960], 
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In  the  computer  program  N  was  allowed  to  range  from  0.2 

to  2.0  i r;  increments  of  0.1.  Values  of  a  less  than  10“^  or 

greater  than  10"1  were  omitted.  The  set  of  alphas  {a..:i<j} 

*  0 

were  computed  for  each  value  of  N.  The  alphas  were  rounded 
off  to  the  first  significant  figure.  Some  of  the  results  of 
these  computations  are  shown  in  Figures  19  through  22.  In 
Figure  19  n  was  set  equal  to  1  and  the  frequency  used  was 
3.15  cps  which  corresponds  to  tnc  predominant  peak  observed 
at  most  distances.  Each  point  P(a,N)  represents  the  number 
of  times  a  particular  value  of  a  occurred  for  a  particular 
value  of  N.  The  actual  numbers  have  been  omitted  here  Tor 
the  sake  of  clarity  but  the  contours  show  the  range  in  values. 
Using  the  same  geometrical  spreading  factor  similar  graphs 
were  made  for  f  *  1.6,  2.0,  2.5,  4.0,  5.0,  6.3,  8,  and  10  cps. 
The  narrow  ridge  of  values  centering  at  a  =  2  x  10"^  was 
founr*  to  shift  slightly  with  frequency.  At  8  to  10  cps,  lar¬ 
ger  values  were  observed  for  N>1  and  a  =  1  x  10"^  while  for 
f  =  1.6  to  2.0  cps  larger  values  of  P  were  fourd  for  N<!  and 
a  =  3  x  1U~^.  The  elongated  ridge  shown  in  Figure  19  was 
somewhat  surprising  since  it  was  thought  that  P  would  have  a 
pronounced  peak  in  the  range  of  N  between  0.7  and  1.0  (Willis 
and  Turpening,  1963). 
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FIGURE  20.  1964  LAKE  SUPERIOR  ATTENUATION  DATA:  a  vs  M.  f  =  3.15  cps 


FIGURE  22.  1964  LAKE  SUPERIOR  ATTENUATION 
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Similar  graphs  were  prepared  for  n  *  1/2,  3/2,  and  2 
at  f  =  3.15  cps.  Two  of  these  are  shown  in  Figures  20  and 
21.  For  the  four  values  of  geometrical  spreading  the  lar¬ 
gest  values  of  P  occurred  at  n  =  3/2.  However,  for  n  =  1, 

P  was  not  significantly  smaller.  A  secondary  peak  in  a 
occurs  at  a  =  10  .  One  would  intuitively  associate  this 

with  the  crustal  attenuation.  However,  these  larger  values 
of  a  are  found  primarily  in  association  with  those  sets  of 
stations  where  the  station  separation  is  small  (A  =  25  km). 

In  those  cases  where  W-  =  W.,  N  would  have  no  significance. 

*  J 

In  order  to  determine  if  this  were  the  reason  for  the  elon¬ 
gated  ridge  at  a  =  2  x  1 0“ 3 ,  all  sets  of  stations  where 
W^/Wj  =  1  were  omitted.  The  results  are  shown  in  Figure  22 
for  n  =  3/2  and  f  =  3.15  cps.  It  can  be  seen  that  there  is 
no  significant  change  in  pattern  from  that  shown  in  Figure 
21.  A  weighting  factor  based  on  the  percentage  of  path  length 
through  the  upper  mantle  compared  to  the  total  path  length 
was  used  for  'the  data  shown  in  Figure  21.  The  significant 
features  remained  unchanged.  Hence,  varying  the  charge  size 
added  an  additional  parameter  that  seriously  interfered  with 
the  interpretation  of  the  attenuation  of  the  seismic  waves 
generated  by  this  series  of  shots.  This  could  have  been  min¬ 
imized  if  one  station  had  been  operated  at  a  small  epicentral 
distance  to  record  all  of  the  shots  at  the  same  point.  In 
this  manner  the  shots  could  have  been  normalized  to  a  stan¬ 
dard  reference  and  corrections  made  accordingly. 
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Early  Rise  Shots 

Figure  23  shows  the  first  arrival  vertical  component  of 
ground  displacement  versus  distance  for  the  Early  Rise  data. 
While  there  is  considerable  scatter,  the  data  show  an  aver¬ 
age  attenuation  proportional  to  the  inverse  3/2  power  of 
the  distance.  On  the  reversed  profile,  the  ECOOE  shots,  the 
1962  North  Carolina  data  and  CHASE  shots  fit  more  closely  to 
an  inverse  third  power  relationship  of  the  distance.  Further¬ 
more,  for  equivalent  sized  shots,  the  shots  fired  in  Lake 
Superior  were  more  efficient  in  coupling  energy  into  seismic 
waves  than  the  shots  fired  in  the  Atlantic  Ocean.  The  sig¬ 
nal  amplitudes  ranged  from  6  to  18  db  or  a  factor  of  2  to  8 
largtr  for  shots  recorded  at  comparable  distances.  A  partial 
explanation  for  this  phenomenon  may  be  that  the  shots  fired 
In  Lake  Superior  (especially  the  Early  Rise  data)  were  fired 
at  an  optimum  depth  for  signal  enhancement  due  to  reinforcement 
by  the  bubble  pulse  (optimum  depth  between  593  and  607  feet). 

Attenuation  measurements  made  in  the  western  United 
States  (Herrin  and  Taggart  [  1  962],  Jordan,  et  al  .  [1965], 

Nuttli  [1963],  Romney  [1959],  and  others)  indicate  a  fall  off 
in  amplitude  proportional  to  the  inverse  cube  of  the  distance. 
Considerable  asymmetry  was  also  disclosed,  especially  in  con¬ 
nection  with  the  GNOME  shot.  Profiles  west  and  northwest 
from  Carlsbad,  New  Mexico  showed  small  amplitude  first  arri¬ 
vals  arriving  several  seconds  late  while  to  the  northeast 
large  first  arrival  amplitudes  could  be  correlated  with  early 
arrivals.  It  is  postulated  that  a  low  velocity  layer  in  the 
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FIGURE  23.  GROUND  DISPLACEMENT  vs  DISTANCE  FOR  PROJECT  EARLY  RISE. 

| 

j 

! 


‘|jw  ^  " 1 " 


upper  mvitle  In  the  western  United  States  could  account  for 
this  anomaly.  The  Lake  Superior  shots  verified  that  there 
was  no  low  velocity  layer  in  the  central  or  eastern  United 
States.  However,  the  shots  fired  in  the  Atlantic  display 
the  inverse  cube  relationship  of  amplitude  fall  off  with 
distance. 

Iso-Particle  Velocity  Maps 

Display  of  seismic  attenuation  data  In  a  reasonable 
manner  has  been  a  difficult  problem  in  the  p'^st.  Presenta¬ 
tions  such  as  shown  In  the  former  section  are  restrictive 
In  that  only  the  peak  signals  are  presented  and  the  frrquen- 
cy  content  is  unknown.  On  the  other  hand  detailed  frequency 
analysis  such  as  shown  In  Figure  18  are  difficult  to  follow 
since  the  spectral  curves  overlap.  The  apparent  scatter  of 
the  data  in  this  form  makes  the  interpretation  very  difficult 
Changes  In  the  spectra  and  rapid  changes  In  amplitude  between 
stations  contribute  to  this  difficulty.  A  method  of  presen¬ 
tation  of  attenuation  data  was  developed  (Willis,  et  al  1964, 
Willis  and  DeNoyer,  1966)  that  permits  presentation  of  the 
attenuation  data  In  three  dimensional  distance,  frequency  and 
particle  velocity  space.  Contour  maps  are  constructed  where 
a  linear  distance  scale  and  a  logarithmic  frequency  scale 
are  used.  The  particle  velocity  values  are  contoured  at  6  db 
Intervals  (a  factor  of  2)  using  a  reference  of  1  x  10"2  in/ 
sec/cycle.  These  amplitudes  correspond  to  true  ground  motion 
All  cf  the  Instrumentation  response  has  been  removed. 
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Figure  24  shows  the  vertical  component  of  the  groun^ 
particle  velocity  for  the  first  three  cycles  of  the  first 
compressional  wave  arrival  recorded  on  the  Early  Rise  shots. 
The  gap  in  data  from  270  to  380  km  is  where  the  profile 
crosses  Lake  Michigan.  The  dots  correspond  to  data  points 
(the  values  have  been  omitted  for  sake  of  clarity).  The 
spectra  of  these  shots  are  strongly  influenced  by  the  bub¬ 
ble  pulse  frequency  and  the  reverberation  of  the  sound 
waves  between  the  lake  bottom  and  the  surface.  The  bubble 
pulse  frequency  is  approximately  2  cps  and  would  be  in  phase 
at  that  water  depth.  The  orqan  pipe  frequencies  for  con¬ 
structive  and  destructive  interference  are: 

Constructi ve  Destructive 

fQ  =  2  cps  2  fQ  =  4  cpr 

3  fQ  =  6  cps  4  fQ  =  8  cps 

5  f Q  =  1 0  cps  6  fQ  =  12  cps 

7  f0  =  14  cps 

A  rather  remarkable  correlation  with  these  frequencies  can 
be  seen  between  the  particle  velocity  high's  and  low’s  shown 
in  this  figure. 

The  major  peaks  along  the  entire  profile  occur  at 
approximately  2  and  5-6  cps.  The  former  coincides  with  the 
bubble  pulse  frequency  and  the  fundamental  organ  pipe  mode 
while  the  latter  coincides  with  the  third  harmonic.  A  high 
can  also  be  correlated  quite  well  with  the  5th  harmonic. 

The  particle  velocity  lows  correlate  quite  well  with  the 
second  and  fourth  harmonic  where  destructive  interference 
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should  occur. 

A  particle  velocity  low  at  all  frequencies  can  be  seer 
near  the  crossover  point  (160  km)  of  the  6.2  and  the  7.7  km/ 
sec  branch  of  the  travel  time  curves.  Unfortunately,  the 
Pn  crossover  point  falls  near  that  section  of  Lake  Michigan 
where  no  d*  ta  are  available.  One  other  predominant  feature 
should  be  discussed.  A  pronounced  particle  velocity  high 
occurs  at  approximately  885  km.  This  same  high  was  observed 
at  a  group  of  four  stations  at  the  same  distance  located 
near  Anna,  Ohio.  The  latter  area  is  150  km  southwest  of  the 
main  profile  and  Is  ar  area  where  small  earthquakes  have 
been  recorded  in  the  past.  The  cause  of  this  large  amplitude 
is  not  known.  It  may  possibly  correlate  with  a  cusp  in  the 
travel  time  curves,  its  bui  dup  is  observable  at  seven  sta¬ 
tions. 

Figure  25  shows  the  1 so-parti cl e-vel ocl ty  map  for  the 
1963  Lake  Superior  South  Line.  The  particle  velocity  high's 
correlated  well  with  the  bubble  pulse  and  organ  pipe  modes 
that  would  cause  reinforcement.  The  predominant  signal  along 
the  entire  profile  is  4  cps  which  correlates  with  the  first 
odd  harmonic  of  the  organ  pipe  mode.  The  sharp  decrease  In 
particle  velocity  amplitude  at  160  to  170  km  correlates  with 
the  data  on  the  Early  Rise  shot  which  occurred  at  the  cross¬ 
over  point  on  the  travel  time  curve  from  the  6.2  to  7.7  branch. 

Other  iso-partlcle-veloclty  maps  of  the  1963  Lake 
Superior  Experiment  were  published  in  a  paper  by  Willis  and 
DeNoyer  [1966],  These  Included  compressional  and  shear  wave 
vertical  component  maps  made  for  data  recorded  on  the 
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Keweenaw  Peninsula  for  shots  fired  in  ihe  eastern  and  wes¬ 
tern  portions  of  the  lake  and  also  for  the  entire  series  of 
shots  recorded  at  Knife  River,  Minnesota.  The  attenuation 
pattern  shown  in  these  maps  is  quite  complex.  Many  of  the 
fluctuations  in  spectral  data  persisted  from  shot  to  shot 
where  the  station  remained  fixed  and  from  station  to  station 
when  the  shot  point  was  fixed.  These  trends  are  difficult 
to  follow  when  examining  the  individual  spectral  curves. 
However,  there  are  many  anomalies  that  can  be  correlated 
at  comparable  distances  on  each  profile.  Likewise  anomalies 
can  be  seen  on  each  line  which  correlate  with  the  shot  and 
are  independent  of  the  epicentral  distance. 

Figures  26  and  27  show  the  particle  velocity  maps  for 
the  first  compressional  wave  arrivals  for  the  1962  North 
Carolina  Experiment.  The  Potter's  Hill  and  Dudley  recording 
stations  were  located  In  line  with  the  shot  profile  and  were 
40  km  *part.  The  similar  major  features  which  correlate  with 
the  same  shot  points  at  these  two  stations  prove  that  con¬ 
ditions  in  the  source  region  strongly  influence  the  charac¬ 
ter  of  the  seismic  signal  recorded  at  distant  stations.  Ver¬ 
tical  and  longitudinal  component  iso-particle-velocity  maps 
for  the  first  P  and  the  maximum  P  arrival  were  made  for  the 
ECOOE  Southern  Profile  150,  250,  350,  400,  450,  and  550/600 
km  stations.  Similar  me  '  r  the  15C,  300,  and  325  km  sites 
on  the  Northern  Profile  were  also  constructed.  Figure  28 
shows  the  vertical  component  particle  velocity  amplitudes 
for  the  150  km  site  on  the  ECOOE  Southern  Profile.  st  of 
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PARTICLE- VELOCITY  MAP  OF  FIRST  P-WAVE  ARRIVAL  FOR  1962  NORTH  CAROLINA  SHOTS 

RECORDEn  AT  DUDLEY,  N.  C. 
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FIGURE  28.  ISO-PARTICLE-VELOCm  MAP  OF  FIRST  P-WAVE  ARRIVAL  FOR  ECOOE  SOUTHERN  PROFILE  SHOTS 

RECORDED  AT  POTTER'S  HILL,  N.  C. 
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these  data  fall  on  the  6.0  branch  of  the  travel  time  curve. 
There  is  not  as  much  relief  or  change  in  signal  level  here 
as  compared  to  the  Early  Rise  or  1963  Lake  Superior  data. 

One  reason  is  that  the  bubble  pulse  frequency  is  below  1.5 
cps  for  all  but  the  two  farthest  shots  and  the  fundamental 
organ  pipe  frequency  is  above  10  cps  for  all  but  the  last 
three  shots.  The  particle  velocity  highs  for  the  last  three 
shots  correlate  with  the  bubble  pulse  frequency  and  the  fun¬ 
damental  mode  organ  pipe  frequency.  All  of  these  maps,  with¬ 
out  exception,  displayed  repeatable  anomalies  that  corre¬ 
lated  with  specific  shot  points. 

Residual  iso-particle-velocity  maps  have  been  prepared 
to  help  isolate  these  anomalies  and  to  determine  the  besl 
values  for  the  absorption  and  geometrical  spreading  losses. 
The  source  spectrum  at  a  standard  reference  distance  was 
computed  by  extrapolating  the  spectral  data  for  the  initial 
P  arrival  computed  at  each  station  to  a  distance  of  1  km 
with  assumed  values  of  a  and  n.  The  average  value  of  the 
particle  velocity  amplitude  was  then  computed.  With  this 
as  the  source  spectrum  the  predicted  spectral  curve  for  each 
shot  was  determined  by  means  of  equation  ?3  with  the  assumed 
values  of  a  and  n.  The  difference  between  the  observed  and 
the  predicted  signal  can  then  be  plotted  and  contoured  simi¬ 
larly  to  the  i so-particl e-vel oci ty  maps.  The  residual  maps 
having  th/*  smallest  overall  values  probably  represent  the 
most  appropriate  values  of  a  and  n  for  the  region.  Examples 
of  these  residual  i  so-part.1  cl  e-vel  oci  ty  maps  may  be  found  in 
the  paper  by  Willis  and  DeNoyer  [1966].  The  correlation 
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between  the  remaining  anomalies,  the  bubble  pul^e  frequency, 
and  the  constructive  and  destructive  interferences  caused  by 
reverberations  of  the  sound  waves  between  lake  bottom  and 
the  water's  surface  Is  quite  plain. 

In  order  to  Investigate  in  more  detail  the  attenuation 
of  the  seismic  waves  using  the  type  of  data  discussed  above 
an  equation  of  the  type 

Al0  "  Aoj  Xin  e"afj  Xl  (25) 

Wo s  employed.  It  accounts  for  geometrical  spreading  (n)  and 
foi  the  absorption  (a)  of  the  seismic  waves.  A^.  represents 
the  particle  velocity  amplitude  for  the  I**1  shot  and  the  j**1 
frequency;  represents  a  constant  assumed  to  be  Independ¬ 
ent  of  distance  but  dependent  on  the  jth  frequency,  repre¬ 
sents  the  distance  from  the  seismometer  to  the  ith  shot;  and, 
fj  Is  the  frequency. 

Initially  we  have  an  equation 


A 


1.1 


1  «  1,2 . M 

j  *  1,2,...,  N 


which  can  be  solveu  for 


eaf  j 


If  we  let 


'oj 


1 


M 


0 

1j 


(26) 


(27) 


(28) 


we  ha»e  arrived  at  equation  (25). 


59 


Hence,  for  each  an  was  calculated  so  that  a  matrix 
was  obtained  whose  rows  (1  ,2, . . . ,j , . . .  ,N)  represent  the  N 
different  frequencies,  and  whose  columns  (1 ,2, . . . ,1 . ,M) 
represent  the  M  different  shot  distances  (in  increasing  or¬ 
der).  An  average  Aqj  was  found,  for  eac!'  frequency,  across 
the  distances  by  summing  the  A.1?  for  a  constant  j  and  divid- 

*  J 

ing  the  sums  by  M,  the  number  of  shots.  There  are  then 
N  AQj's  which  are  in  a  respect  independent  of  distance. 

This  represents  an  average  spectrum  (Aqj). 

Using  equation  (25)  and  the  computed  AQj  ,  new  A^ j's 
designated  Afj)  are  computed,  which  in  general  will  be  dif¬ 
ferent  from  the  original  A ^ j ' s . 

The  A^'s  and  Afj's  are  converted  from  particle-velocity 
amplitudes  to  the  equivalent  amplitudes  in  decibels  by  using 
the  following  relationship: 

An  (db) 

A.  (in/sec)  =  20  log10  -  (29) 

1J  lu  0.01 

A  matrix  cf  residualr  S..)  are  formed  where 

*  J 


(30) 


The  S..'s  are  plotted  on  a  graph  where  the  abscissa 

*  J 

corresponds  to  the  distance  and  the  ordinate  is  the  frequency. 

The  summation  of  the  residuals  vs  the  aL^orption  co¬ 
efficient  for  values  of  n,  the  geometrical  spreading  factor, 
ranging  from  0  to  3  are  shown  in  Figure  29  for  the  Lake 
Superior  1963  South  Line.  The  best  fit  for  n  and  a  would 
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SUMMATION  OF  RESIDUALS  vs  a  FOR  1R63  LAKE  SUPERIOR  SOUTH  LINE 
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correspond  to  the  residual  minimum.  For  this  profile  a 
value  of  n  =  2.5  and  a  =  7  x  10"^  yields  the  smallest  re¬ 
sidual.  A  value  of  n  =  2.0  would  correspond  closely  to  the 
geometrical  spreading  due  to  he'd  wave  propagation.  For 
n  =  2.0,  a  =  1  x  10“^  gives  the  best  fit.  This  would  cor¬ 
respond  to  a  Q  of  approximately  490. 

Q  refers  to  the  seismic  anelasticity  of  a  medium  and 
2ttE 

can  be  defined  as  where  E  is  the  energy  in  the  seismic 

wave.  An  equivalent  expression  given  by  Haskell  r  057]  de_ 
fines  Q“1  as  26  x  velocity  which  in  the  notation  of  this 
paper  would  be  — — oc  i  ty  Spectl  um  0f  a  seismic 

signal  is  a  function  of  the  source,  the  seismograph,  and 
the  transmission  path.  If  the  source  function  is  known  or 
is  a  constant  and  the  response  of  the  seismograph  is  known 
any  systematic  variations  in  the  seismic  signal  can  be 
correlated  with  transmission  path  affects.  As  mentioned 
earlier  these  affects  are  due  to  geometrical  spreading  and 
to  anelastic  absorption  or  scattering.  The  techniques  des¬ 
cribed  above  were  an  attempt  to  determine  the  best  fit  for 
the  geometrical  spreading  factor  and  the  absorption.  De¬ 
pending  on  wave  type,  the  geometrical  loss  correction  fac¬ 
tor  should  be  easily  removed  from  the  attenuation  data  (n  - 
0.5  for  surface  waves,  n  =  1  for  body  waves  and  n  =  2  for 
head  waves).  However,  it  was  found  that  for  first  compres- 
sional  wave  arrivals,  values  of  n  other  than  1  or  2  yielded 
smaller  residuals  in  many  instances.  This  variability  in 
n  strongly  influences  the  values  one  would  obtain  for  the 
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absorption  coefficient.  The  Early  Rise  shots  were  the  only 
ones  where  the  source  spectrum  remained  constant.  Hence, 
caution  should  be  used  In  comparing  the  values  of  Q  from 
the  different  profiles. 

The  attenuation  data  for  the  1963  Lake  Superior  South 
Line  are  somewhat  anomalous  compared  to  the  other  profiles 
described  In  this  report.  However,  If  only  the  data  from 
the  predominant  frequencies  are  used  (4  to  8  cps)  the  disa¬ 
greement  Is  reduced  significantly.  At  4  cps  the  best  fit 
for  this  profile  Is  given  by  a  value  of  n  =  1.0,  a  =  4  x  10-3 
(Q  »  123),  and  n  *  1.5,  o  «=  3  x  10"3  (Q  -  164).  On  the 
1963  Lake  Superior  eait  line  and  Knife  River  profiles  we 
have  values  of  n  *  1.5,  o  a  3  x  10“3  (Q  *  158)  and  n  =  1.5, 
a  =  5  x  10'4  (Q  *  894)  respectively.  These  minima  agree 
closely  with  the  minima  obtained  using  all  of  the  frequency 
data.  It  would  thus  appear  that  local  seismic  background 
noise  adversely  affected  portions  of  the  spectrum  where  the 
signal -to-nolse  level  was  low  on  the  south  line  introducing 
a  bias  Into  the  attenuation  studies. 

The  Early  Rise  data  were  analyzed  In  a  similar  manner. 
Spectral  analyses  were  limited  to  recordings  made  at  dis¬ 
tances  less  than  1068  km  due  to  the  poor  s ignal -to-nolse 
levels  at  greater  distances.  The  vertical  component  of  the 
first  compressi onal  wave  arrival  over  this  distance  range 
had  a  minimum  residual  at  n  ■  1  and  a  ■  2  x  10"4  (Q  =  1895). 
For  the  second  branch  of  the  travel-time  curve  corresponding 
to  the  7.7  km/sec  layer  the  minimum  residual  occurred  at  a 
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value  of  n  =  1.5  and  a  =  2  x  IO"3  (Q  =  204). 

Figure  30  shows  the  sum  of  the  residuals  versus  a  for 
n  ranging  from  0.5  to  3  for  th:.  ECOOE  southern  profile  150 
km  site.  The  best  fit  is  for  n  =  0.5  and  a  =  6  x  10"^  (Q  = 
873).  For  n  =  1.0,  a  value  of  a  =  3  x  IO"4  yields  the 
smallest  residuals.  It  can  be  seen  in  this  figure  that 
values  of  N  greater  than  1.5  can  be  discarded.  However, 
even  though  the  value  of  n  =  0.5  gives  a  smaller  residual 
than  n  =  1.0,  the  difference  is  not  large  so  there  is  some 
uncertainty  which  value  of  n  is  most  appropriate. 

Similar  graphs  were  prepared  for  the  other  recording 
stations  of  the  ECOOE  series.  The  average  values  of  a  for 
the  minimum  residuals  are  summarized  in  Table  III. 
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The  Lake  Superior  data  showed  a  degree  of  consistency 
with  minimum  residuals  corresponding  to  values  of  n  =  1.5 
and  Q  varying  from  158,  164,  204,  and  894  for  crustal  paths. 
The  latter  (Knife  River  profile)  also  included  waves  which 
penetrated  the  upper  mantle.  The  Early  Rise  profile  which 
included  long  transmission  paths  through  the  mantle  yielded 
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a  value  of  Q  equal  to  1895.  The  east  coast  data  showed  lar¬ 
ger  crustal  Q  values  and  smaller  upper  mantle  Q  values  than 
the  Lake  Superior  data.  Intuitively  one  would  anticipate 
more  reliable  values  of  Q  for  the  Lake  Superior  data  since 
the  east  coast  data  are  Influenced  by  the  Appalachian  root 
system  and  the  transition  from  oceanic  to  continental  crusts 
and  mantles. 

The  range  In  values  of  Q  reported  here  Is  not  inconsis¬ 
tent  with  published  values.  Anderson,  et  al  [1965]  esti¬ 
mates  that  tne  upper  mantle  Q  Is  approximately  1500  for  com- 
presslonal  waves.  Sumner  [1967]  reports  values  of  Q  greater 
than  1000  for  compresslonal  waves  In  the  Andes.  O'Brien 
[1967]  and  Mueller  [1968]  give  values  of  crustal  Q's  of  830 
±30  and  670  ±100  respectively.  Teng  [1968]  reports  values 
of  Q  of  450  in  the  upper  mantle  with  a  low  Q  zone  of  60  at 
a  depth  of  50  km  below  the  upper  mantle  boundary.  In  the 
western  United  States  Long  [1968]  has  found  values  of  Q  at 
distances  between  100  and  600  km  of  116  for  the  SHOAL  event 
and  169  for  the  GNOME  event.  Other  Investigators  (Anderson 
and  O'Connell  [1967],  Anderson  [1967],  Anderson  [1966],  and 
Walsh  [1968])  have  studied  shear  and  surface  wave  Q  values 
In  the  mantle  as  a  function  of  depth  and  temperature.  Their 
studies,  however,  are  restricted  primarily  to  long  period 
data.  Average  values  of  Q  for  surface  waves  In  the  upper 
400  km  of  the  mantle  range  between  88  and  135.  Shear  Q 
values  average  about  200  for  the  upper  600  km  of  the  mantle 
and  about  600  for  the  entire  mantle. 
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Pakiser  [1965]  has  shown  that  the  Rocky  Mountains  di¬ 
vides  the  crust  and  upper  mantle  of  the  United  States  into 
two  superprovinces.  The  eastern  United  States  has  a  thick¬ 
er  average  crust  with  larger  crustal  and  Pn  velocities. 

This  region  also  has  larger  magnetic  anomalies  indicating 
an  abundance  of  magnetic  materials.  The  higher  velocities 
in  the  lower  crust  could  be  explained  by  lateral  variations 
of  ec.ogite  to  amphibolite.  The  average  upper  silicic 
crust  is  about  the  same  thickness  (19  km)  in  the  eastern 
and  western  United  States  while  the  average  lower  mafic 
crust  is  thicker  in  the  eastern  United  States. 

The  higher  crustal  and  mantle  P  velocities  of  the 
eastern  United  States  correlate  with  the  larger  values  of 
Q  found  in  this  same  area.  The  smaller  upper  mantle  Q 
values  observed  in  the  western  United  States  are  probably 
due  to  the  low  velocity  layer  in  the  upper  mantle  although 
it  may  also  be  due  in  part  to  lateral  variations  of  eclogite- 
peridotite  in  the  upper  mantle.  Since  high  temperatures 
cause  the  attenuation  of  seismic  waves  to  increase  rapidly 
(Anderson  [1967.1),  the  higher  heat  flow  measurements  in  the 
western  United  States  may  also  have  some  correlation  to  the 
lower  Q  values. 

In  Figure  31  the  summation  of  residuals  versus  n  is 
presented  where  the  absorption  coefficient,  a,  was  set  equal 
tc  zero.  All  of  the  attenuation  of  the  seismic  waves  thus 
is  accounted  for  by  the  geometrical  spreading  factor.  The 
minima  vary  from  n  =  1  (ECOOE)  to  n  =  3.5  (Lake  Superior 
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FIGURE  30 


.  SUMMATION  OF  RESIDUALS  vs  a  FOR  ECOOE  SOUTHERN  PROFILE 
150  km  SITE. 


FIGURE  31.  SUMMATION  OF  RESIDUALS  vs  n. 
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South  Line)  and  average  approximately  1.5.  Hence,  geometri¬ 
cal  spreading  by  head  waves  (n  =  2.0)  would  cause  too  great 
a  loss  in  eisrnic  energy  to  account  for  the  observed  loss 
at  most  of  these  sites.  The  head  wave  amplitudes  were  ei¬ 
ther  too  small  to  be  observed  or  perhaps  the  upper  mantle 
boundary  is  too  uneven  to  support  their  transmission. 


Additional  Studies 


The  Investigation  of  the  atter*-  tlon  of  seismic  waves 
discussed  In  this  report  did  not  take  Into  consideration 
the  effect  of  the  changing  anale  of  Incidence  of  the  ceismic 
wave  as  a  function  of  distance  from  the  sou/ce.  Three 
component  data  were  obtained  for  all  of  the  shots  reported 
here.  Particle  motion  diagrams  have  been  prepared  for  the 
first  compresslonal  wave  arrivals.  A  linear  least-square 
regression  analysis  technique  was  Investigated  for  estimat¬ 
ing  the  parameters  affecting  the  attenuation  of  seismic 
waves.  A  negative  absorption  coefficient  for  the  longitudi¬ 
nal  component  was  disclosed  by  the  least  square  approach 
which  Implies  that  the  vertical  and  longitudinal  component 
particle  velocities  should  be  combined  vectorlally  and  that 
the  transformed  data  should  yield  better  estimates  of  the 
propagation  losses. 

The  problem  of  shot  point  and  recording  station  anoma¬ 
lies  should  be  Investigated  further  to  determine  the  best 
approach  for  removing  these  affects  in  computing  the  absorp¬ 
tion  coefficients.  The  former  anomalies  are  not  too  criti¬ 
cal  where  the  same  shot  point  Is  used  (e.g.,  the  Early  Rise 
shots^.  However,  where  moving  shot  points  and  fixed  record¬ 
ing  stations  are  used  {e.g.,  1963  Lake  Superior  Experiment 
and  ECOOE),  the  shot  point  anomalies  can  frequently  mask  or 
seriously  confuse  the  normal  attenuation  losses. 

Another  topic  that  needs  further  Investigation  Is  the 
particularly  good  coupling  of  seismic  energy  found  In  the 
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Lake  Superior  area.  Current  studies  are  being  conducted 
that  nu,y  give  some  insight  into  this  subject.  A  series  of 
recordings  were  made  in  Lake  Superior  on  the  lake  bottom 
near  shot  point  51  during  the  summer  months  of  1967.  These 
recordings  were  made  w*th  modified  ocean  bottom  seismographs 
in  an  effort  to  determine  if  there  was  any  reciprocity  in 
the  propagation  of  seismic  waves  from  teleseismic  events. 
These  recordings  are  currently  being  analyzed. 

Cone! usions 

1.  Seismic  refraction  studies  disclosed  that  the 
earth's  crust  varies  in  thickness  from  28.1  km  at  the  con¬ 
tinental  margin  in  North  Carolina  to  50.4  km  near  shot  point 
51  in  Lake  Superior.  A  single-layer  crust  with  a  congres¬ 
sional  velocity  of  6.0  km/sec  was  found  in  North  Carolina 
while  a  two-layered  crust  with  compressi onal  velocities  of 
6.2  and  7.7  km/sec  was  found  in  the  Lake  Superior  area. 
Larger  travel  time  residuals  for  recording  stations  located 

i  the  Appalachian  Mountains  are  indicative  of  a  mountain 
root  system  similar  to  the  Rocky  Mountains.  No  indication 
of  a  low  velocity  layer  was  found  in  the  upper  mantle  in  the 
eastern  United  States  that  would  correlate  with  such  a  fea¬ 
ture  reported  in  the  western  United  States. 

2.  The  Lake  Superior  area  was  found  to  be  more  effi¬ 
cient  in  the  coupling  of  seismic  energy  for  underwater  shots 
than  the  area  in  the  Atlantic  Ocean  off  the  coast  of  North 
Carolina.  Signals  ranged  up  to  a  factor  o'  eight  larger. 
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3.  A  new  method  of  displaying  seismic  attenuation  and 
spectral  data  was  developed  that  permits  the  presentation 

of  the  data  In  three-dimensional  space — distance,  frequency, 
and  ground  particle  velocity  amplitude.  This  technique  al¬ 
lows  the  presentation  of  large  *  ^entitles  of  data  in  a  con¬ 
cise  manner  that  assists  In  the  interpretation  of  the  data. 

4.  Meaningful  and  consistent  anomalies  in  the  spec¬ 
trum  of  seismic  waves  generated  by  underwater  shots  were 
found  that  could  be  correlated  with  variations  In  the  source 
region  and  In  propagation  effects. 

5.  The  maximum  amplitude  of  the  first  compressl onal 
wave  arrivals  attenuate  like  the  Inverse  3/2  powe-  of  the 
dlstanc*  out  to  approximately  1100  km  from  the  Lake  Superior 
area.  Conversely  an  Inverse  cube  law  applies  to  shots  fired 
in  th2  Atlantic. 

6.  The  detailed  spectral  studies  of  the  Lake  Superior 

shots  disclosed  that  the  crustal  attenuation  data  could  be 
fit  using  a  geometrical  spreading  factor  of  n  *  3/2  and  a  s 
2-3  x  10“3.  Mantle  compressl onal  waves  fit  a  geometrical 
spreading  function  of  n  -  1  and  a  *  2  x  10"4.  These  give 
average  upper  crustal  Q  values  of  approximately  160,  a  low¬ 
er  crustal  Q  of  approximately  200  and  an  upper  mantle  Q  of 
approximately  1900.  The  Atlantic  Ocean  shots  had  values  of 
n  *  12  and  o  *  6  10"4  to  4  x  10“^.  The  corresponding 

values  of  Q  would  equal  873  and  960.  The  former  includes 
transmission  paths  predominantly  within  the  crust  while  the 


71 


latter  has  transmi  ssi on  paths  predominantly  in  the  mantle. 
Hence,  the  east  coast  data  shows  the  same  order  of  magnitude 
values  for  Q  in  the  crust  and  upper  mantle. 

7.  Attenuation  losses  in  the  propagation  of  seismic 
waves  did  not  support  the  theory  of  head  wave  propagation 
which  would  require  a  decrease  in  amplitude  approximately 
proportional  to  the  inverse  square  of  the  distance. 

8.  The  variation  in  values  for  the  geometrical  spread¬ 
ing  factor  and  absorption  coefficient  caused  by  the  broad 
minimum  residuals  indicate  that  scattering  effect-,  station 
anomalies  and  source  variations  play  an  importart  role  in 
the  apparent  attenuation  of  seismic  waves  and  that  methods 
of  determining  the  influence  of  these  factors  will  have  to 
be  developed  before  more  accurate  determi nations  of  Q  can 

be  made  over  the  frequency  range  of  .5  to  25  cps. 
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APPENDIX  A 

Underwater  Explosion  Experiments 

The  series  of  underwater  shot  programs  listed  in  Table 
I  were  participated  in  by  a  number  of  organizations  that 
banded  together  in  the  loosely  knit  group  called  the  North 
American  Explosion  Seismology  Group.  Participation  of  the 
members  varied  with  each  series  of  shots.  The  participants 
are  listed  below. 


Boston  College 

Carnegie  Institution  of  Washington 
Dominion  Observatory,  Canada 
Georgia  Institute  of  Technology 
Geotech  Division,  Teledyne  Industries 
Pennsylvania  State  University 
Princeton  University 
Stanford  Research  Institute 
Southwest  Center  for  Advanced  Studies 
University  of  Alberta 
University  of  Manitoba 
University  of  Michigan 
University  of  Minnesota 
University  of  Toronto 
University  of  Western  Ontario 
University  of  Wisconsin 
U.  S.  Geological  Survey 
Vela  Sei smol ogl cal  Center,  AFTAC 

The  1962  North  Carolina  Experiment  i.as  conducted  during 
the  months  of  July  and  August.  The  program  consisted  of  a 
series  of  39  one-ton  shots  fired  on  the  ocean  bottom  along  a 
profile  normal  to  the  coast  near  Cape  Fear,  North  Carolina. 


The  shots  were  fired  on  and  off  the  continental  shelf.  De¬ 
tails  of  this  experiment  may  be  found  in  a  paper  by  Meyer, 
et  al .  [1965]. 

The  1963  Lake  Superior  Experiment  consisted  of  a  series 
of  78  one-ton  shots  fired  along  a  line  extending  the  length 
of  Lake  Superior.  All  shots  were  fired  on  the  lake  bottom. 
Details  of  this  experiment  have  been  described  by  Steinhart 
[1964]. 

The  October  1964  Lake  S  'perior  Experiment  consisted  of 
a  series  of  1 ,  2,  3,  6,  and  10  tons  of  Nitramon  WW-EL  fired 
on  the  bottom  at  shot  points  45-49  and  73  of  the  1963  series. 
A  total  of  10  shots  were  fired.  One  shot  of  each  size  was 
fired  at  shot  point  45. 

The  East  Coast  On-Shore  Off-Shore  Experiment  (EC00E) 
consisted  of  a  series  of  one-,  five-,  and  ten-ton  shots 
fired  along  four  lines,  two  off  the  coast  of  North  Carolina 
and  two  off  the  coast  of  Virginia.  Details  of  this  experi¬ 
ment  may  be  found  in  a  report  by  Hales,  et  al  [1967]. 

The  last  large  scale  cooperative  experiment  consisted 
of  the  38  five-ton  shots  of  the  EARLY  RISE  series.  These 
were  fired  in  Lake  Superior  on  the  lake  bottom  near  shot 
point  51  of  the  1963  series.  Details  of  this  program  are 
in  a  report  by  Warren,  et  al  [1967]. 

The  CHASE  (cut  holes  and  sjnk  em)  series  consisted  of 
three  shiploads  of  surplus  explosives  that  were  disposed  of 
by  the  U.  S.  Navy.  The  explosive-laden  ships,  primed  with 
pressure  detonators,  were  scuttled  in  deep  water  off  the 
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Atlantic  Coast.  Pertinent  details  may  be  found  in  Table  IV. 
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APPENDIX  B 

Instrumentation  Characteristics 


Matched  three  component  seismometers  of  the  Hall-Sears 
Model  HS-10  2  cycle  type  were  used  for  most  of  the  data  pre 
sented  in  this  report.  High  Impedance  colls  (200,000  ohms) 
were  used  for  high  gain  purposes.  The  transistorized  seis¬ 
mic  amplifiers  used  were  constructed  by  the  Geophysics  Lab- 
orabory  and  had  a  maximum  gain  of  114  db.  The  bandwidth  of 
these  Instruments  was  essentially  flat  from  0.2  to  500  cps 
and  had  an  Input  noise  with  a  5  megohm  source  of  less  than 
1  yV.  A  variety  of  Instrumentation  type  FM  tape  recorders 
were  used  that  had  an  average  signal -to-nol se  ratio  of  40 
db.  The  standard  practice  used  In  the  field  was  to  record 
the  output  of  each  seismometer  at  two  or  three  levels 
separated  by  12  or  18  db  to  provide  wider  dynamic  range.  A 
more  detailed  description  of  this  equipment  may  be  found  In 
Willis  and  Jackson  [1966],  Willis  and  Johnson  [1959],  and 
Willis  [1966]. 
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APPENDIX  C 

Spectral  Analyses  Techniques 


The  frequency  analyses  data  presented  in  this  paper 
were  obtained  by  bandpass  RLC  filtering.  Different  spec¬ 
tral  analyses  techniques  were  investigated.  These  include 

Bandpass  RLC  filtering 

Fourier  transform  analysis  employing  an 
electronic  analog  computer 

Digital  analyses  using  an  IBM  7090  em¬ 
ploying  Fourier  transform  methods 

Optical  diffraction  of  variable-density 
film  records. 

A  comparison  of  the  first  three  techniques  may  be  found 
in  a  paper  by  Mclvor  [1964].  The  optical  analysis  technique 
is  described  in  papers  by  Jackson  [1964]  and  Jackson  [1965]. 
The  results  obtained  from  these  different  approaches  are  in 
good  agreement  if  the  data  are  normalized  to  the  same  band¬ 
width.  The  digital,  analog  and  optical  computer  technique 
provide  finer  frequency  resolution  than  the  bandpass  filter¬ 
ing  method.  However,  for  low  frequency  seismic  signals  (be¬ 
low  10  cps)  this  fine  frequency  resolution  usually  is  not 
required. 

For  the  data  reported  in  this  report  the  spectral  in¬ 
formation  was  performed  by  playing  back  the  demodulated  FM 
magnetic  tape  signals  through  a  series  of  electronic  poly¬ 
filters,  each  filter  of  which  was  spaced  at  one- thi rd-octave 
intervals.  These  filters  are  composed  of  three  transis¬ 
torized  cascaded  RC  networks  with  active  feedback  elements. 
The  outputs  are  passed  through  DC  line  amplifiers  to  bifilar 
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galvanometers  with  a  natural  frequency  of  500  cps.  The 
galvanometer  deflections  are  recorded  on  a  conventional 
photographic  recorder.  The  time  window  of  these  filters 
Is  approximately  Gaussian.  The  use  of  the  DC  line  ampli¬ 
fier  permits  amplification  of  the  signal  thus  facilitating 
the  Interpreters  measurement  of  the  trace  deflection  (pro¬ 
viding  of  course,  that  the  signal  Is  not  In  the  seismic 
background  noise  or  the  noise  of  the  recording  and  play¬ 
back  systems).  This  permits  more  accurate  measurements  and 
In  many  cases  allows  the  Interpreter  to  measure  signals 
that  would  otherwise  be  omitted  because  of  low  signal  levels. 
One  can  thus  examine  teleseismlc  signals  that  contain  energy 
well  below  the  natural  frequency  of  the  seismometer.  This 
Is  Illustrated  in  Table  V. 

Shown  In  Table  V  are  the  signal  levels  (*)  for  the  first 
compressi onal  wave  arrival  of  a  teleseismlc  event  as  meas¬ 
ured  on  each  one-third-octave  seismogram  over  the  frequency 
range  from  0.315  to  2.5  cps.  At  0.4  cps  the  signal  level 
read  on  the  third-octave  seismogram  was  -10  db  (peak-to-peak 
tract  deflection  of  0.632  inches).  If  no  playback  gain  were 
used,  this  signal  level  would  have  been  -30  db  (peak-to-peak 
trace  deflection  of  0.063  Inches).  The  signal  level  was 
10  db  (factor  of  3.2)  above  the  background  noise.  Hence, 
it  is  evident  that  It  is  much  more  accurate  to  read  a  dis¬ 
placement  of  0.632  inches  with  an  average  background  noise 
of  0.2  Inches  than  0.063  Inches  with  a  background  of  0.02 
Inches.  A  signal  of  this  amplitude  would  often  be  overlooked 
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especially  if  one  were  measuring  the  signal  level  on  the 
next  lower  frequency  seismogram  (0.315  cps). 

Column  (a)  shows  the  corrected  seismometer  signal  out¬ 
put  which  includes  corrections  for  the  record  gain  used 
in  the  field  recording  and  the  playback  gain  used  in  the 
laboratory.  Column  ( b )  takes  into  account  corrections  for 
the  seismic  preamplifier  response  and  the  seismometer  re¬ 
sponse.  The  signal  levels  shown  here  would  be  for  an  am¬ 
plifier  and  a  seismometer  whose  response  would  be  flat  be¬ 
tween  0.315  and  2.5  cps.  Column  (a)  shows  a  peak  signal 
level  at  0.8  cps.  At  0.5  cps  the  signal  level  is  11  db  low¬ 
er.  However,  from  column  (b)  it  can  be  seen  that  after 
corrections  for  instrumentation  response  are  made,  the  dif¬ 
ference  in  signal  level  between  0.8  and  0.5  cps  is  only  1  db. 
Hence,  the  true  energy  has  a  broader  spectral  width  than 
is  apparent  on  a  broadband  seismogram.  The  peak  signal  is 
below  the  seismometer's  resonant  frequency  and  hence  any 
variation  in  frequency  of  this  signal  from  one  recording 
site  to  the  next  would  be  magnified  or  exaggerated  on  a 
broadband  or  regular  seismogram  (for  example,  see  Willis 
1965).  This  variation  would  not  be  as  evident  or  might  be 
completely  masked  on  the  latter  type  record. 
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Table  V.  Spectral  Corrections 
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eismometer  response  corrected  to  a  flat  output  of  1  Vrms  at  a  constant  velocity 
10-2  inches  per  second. 
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